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pH , Abstract. A general method is discussed to obtain Markovian master equations which 

O |' describe the interaction with the environment in a microscopic and non-perturbative fashion. It 

is based on combining time-dependent scattering theory with the concept of continuous quantum 
measurements. The applications to the case of a Brownian point particle and to the case of a 
' complex molecule, both in the presence of a gaseous environment, are outlined. 
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1 1. Introduction 

Quite a number of contributions presented at the DICE2006 workshop were focused on possible 
extensions of standard quantum theory — be it to overcome the incompatibility of quantum 
theory and general relativity, or to understand the measurement process and our perception 
^2 , of the world as classical. When discussing these modifications of quantum theory, which often 

imply a violation of unitarity, it is important to keep in mind that their observability may 
be severely limited by environmental decoherence pQ. The latter occurs naturally within the 
framework of quantum theory if one accounts for the fact that all practically relevant quantum 
systems are in contact with some uncontrollable environment. Before interpreting the effects 
and judging the relevance of non-standard extensions it is therefore of great importance to have 
a full quantitative understanding of environmental decoherence. 

In the present contribution I would like to discuss a new method to derive Markovian 
master equations for open quantum systems, which take into account the interaction with the 
environment in a microscopic and non-perturbative fashion [2j [3] . Such dynamic equations are a 
prerequisite for any faithful quantitative description of open quantum systems. In particular, it 
is well-known that phenomenological or perturbative master equations may fail by many orders 
of magnitude when applied to situations where the decoherence time is much shorter than the 
typical time scale of dissipation, even though they describe dissipation phenomena rather well, 
see e.g. J4|. 

In the standard approaches one starts out with an approximate total Hamiltonian for the 
system plus environment [5j. Then, various approximations are employed in the course of 
the derivation, one of which is often the Markov assumption. It implies that environmental 
correlations disperse fast, so that on a coarse-grained timescale the temporal change of the system 
state depends on the present state of the system, but not on its history. The starting point of 
the method presented here differs considerably from the standard methods. It is not based on 
an approximate "total" Hamiltonian of system plus environment. Rather, the environmental 




Figure 1. (a) It is assumed that the system interacts at most with one environmental (quasi-) 
particle at a time (no three-body collisions) and that they disperse their correlation with the 
system before scattering again. This is tantamount to the Markov assumption, (b) In order to 
consistently incorporate the state-dependence of the collision rate into the dynamic description, 
we imagine that the system is monitored continuously by a transit detector, which tells with a 
temporal resolution At whether a particle is going to scatter off the system, or not. 

coupling is described in an operational sense, by using scattering theory, which permits a non- 
perturbative description. In particular, this admits to take the Markov assumption as a premise, 
rather than having to introduce it 'by hand' in the course of the calculation. 

The use and strength of the method will be exemplified by discussing the dynamics of two 
different quantum systems in the presence of an ideal, thermalized gas — a trapped molecule 
and a Brownian point particle. The corresponding master equations will be given in terms of 
the exact scattering amplitudes describing the interaction of the gas particle with the quantum 
system. 

2. The monitoring approach 

The present approach to the derivation of master equations is motivated by the observation that 
environmental decoherence, and Markovian open quantum dynamics in general, can be viewed as 
due to the information transfer from the system to the environment, see e.g. [3]. In accordance 
with this, we will picture the environment as consisting of (quasi-)particles, which serve to 
monitor the system continuously in terms of microscopic probes scattering off the system at 
random times. This point of view will be applicable whenever one can describe the interaction 
with the environment in terms of individual interaction events, and, of course, it suggests a 
formulation in terms of scattering theory. The Markov assumption is then easily incorporated 
from the outset by disregarding the change of the environmental state after each collision. 

When setting up a dynamic equation, one would like to write the change of the system 
state in time as the infinitesimal rate of scatterings multiplied by its effect, i.e., by the state 
transformation due to an individual collision. However, not only the transformed state depends 
on the original system state, in general, but also the collision rate, so that this naive ansatz would 
yield a nonlinear equation. In order to account properly for the state dependence of the collision 
rate we will make use of the concept of continuous and generalized quantum measurements 
[H [8]. Specifically, and just for the sake of consistent bookkeeping of the probabilities, we 
shall assume that the system is surrounded by a hypothetical, minimally invasive transit detector, 



which tells at any instant, and with a time resolution At, whether a probe particle has passed 
by and is going to scatter off the system, see Fig. [TJ 



2. 1 . The effect of continuous monitoring 

Let us denote by C the event that the transit detector triggers in the interval (t; t + At) 
indicating that a collision is imminent. In general, the rate of collisions is then described by a 
positive operator V acting in the the system-probe Hilbert space. Given the uncorrelated state 
£tot = P ® Pe h determines the probability of a collision to occur in a small time interval At, 

Prob(C in At\p®p E ) = Attr (T [p <g> p E ]) . (1) 

Here p E is the stationary reduced single particle state of the environment. The microscopic 
definition of V will in general involve the current density operator of the relative motion and a 
total scattering cross section, see below. 

It is now important to note is that the information that a collision is going to take place 
changes our knowledge about the state, and therefore its description. Most generally, it is 
transformed by a generalized measurement transformation, i.e., a trace decreasing, completely 
positive map [7j- At the same time, we have to take into consideration that the detector is 
not real, but is introduced here only for enabling us to account for the state dependence of the 
collision probability. It is therefore reasonable to take the detection process as efficient and 
minimally-invasive, so that it is described by a single operator which introduces no reversible 
back-action. The consistency requirement for measurement transformations then implies that 
after a (hypothetical) detector click, but prior to scattering, the normalized system-probe state 
must have the form 

rV2 , .ri/2 

This nonlinear transformation reflects our improved knowledge about the incoming two-particle 
wave packet, and it may be viewed as enhancing those parts which are heading towards 
a collision. Also the absence of a detection event during At constitutes a measurement 
which changes the state. The corresponding probability must satisfy Prob [C in At\g^ t) = 
1 — Prob (C in At\g^ Q ^), and the state conditioned on a null-event is given by 

Wotl I 1- A(tr(r ftot ) 

In fact, more general nonlinear transformations are conceivable, but this one is distinguished by 
the fact that it introduces no further operators. 

The effect of a single collision is described by the two-particle S-matrix S. The unconditioned 
system-probe state after time At can now be formed by taking into account that the detection 
outcomes are not really available. The infinitesimally evolved state is then given by the sum 
of the colliding state transformed by the S-matrix and the untransformed non-colliding one, 
weighted with their respective probabilities, 



g[ ot (At) = Prob (C in At\g tot ) SM (g tot \C) S+ + Prob (C in At\g tot ) M {g tot \C) 

= srVW^stA* + etot - r^ tot r 1/2 At. 



(4) 



2.2. Constructing a dynamic equation 

In order to obtain a differential equation for the time evolution it is now convenient to split off 
the nontrivial part of the two-particle S-matrix, thus defining 



T = i(I-S). (5) 

This operator is proportional to the T-matrix of scattering theory (only) on the energy shell, 
and one finds that the unitarity of S implies that it satisfies i (T — T^) = — T^T. Using this 
relation the differential quotient can be written as 



g/(A f } = JY^qY 1 ' 2 ^ - -VjY^gY 1 ' 2 - ir x /V 1/2 TtT + + V, Y^gY 1 / 2 
At y 2 u 2 e 2 L 5 ^ 



(6) 



From here it is easy to obtain a closed differential equation for p. We trace out the environment 
with q = p® p E , take the limit of continuous monitoring At — > 0, to arrive at [3] 

^t P = ^P] + ^E([T + T\V 1 / 2 [p®p E ]Y 1 / 2 ])+Tr E (jY'/ 2 [p®p E ]Y 1 / 2 ^ 

-I-Tve ( rV2 T t Tr i/2 [p g, Pe] ) _ J tv s ( [p 9 Pe] Y^VJY 1 ' 2 ) . (7) 



The Markov approximation enters here by assuming the factorization p ® p E to be valid for all 
times. Note also that the generator H of the free system evolution was added, thus switching from 
the interaction picture to the Schrodinger picture. The collision rate with its state dependence 
is incorporated by the operators l~ 1//2 , while the operators T describe the individual microscopic 
interaction process without approximation. 

The discussion was very general, so far. However, to obtain concrete master equations 
system and environment have to be specified, along with the operators Y and S describing 
their interaction. In the following applications, we will assume the environment to be an ideal 
Maxwell gas. Its single particle state 

A ?h... J »p 2 



is characterized by the inverse temperature j3, the normalization volume £1, and the thermal de 
Broglie wave length = ^J2irh 2 f3 /m. 



3. Collisional decoherence of a discrete object 

As a first concrete implementation of the monitoring approach, let us see how an immobile 
object with discrete internal states, such as the electronic, configuration, and rotation dynamics 
of a trapped molecule, gets affected by an environment of ideal gas particles. For example, 
such a master equation will describe how a possible superposition of configuration eigenstates 
of macromolecules decoheres into a mixture of robust configuration states, see Fig. [21 

The interaction between system and gas will be described in terms of the exact scattering 
amplitudes f aao (p,p ) determined by the interaction potential. They describe the in general 
inelastic coupling between the internal energy eigenstates, which form a discrete basis of the 
system Hilbert space and are called channels in the language of scattering theory. The notation 
| a) will be used in the following to indicate the system eigenstates of energy E a . In this channel 
basis, p a p = (a\p\(3), the equation of motion ([7|) takes on the form of a general discrete master 
equation of Lindblad type, 

d tPaP = -[u; a pp^+^p ao p M^f -^Y,P^M^ a -^p a p M^°. (9) 

aoPo "07 A)7 




Figure 2. Superpositions of 
different configuration states of 
a molecule will decohere into a 
mixture of configuration states, as 
a result of the interaction with a 
gaseous environment. The master 
equation ([9]) describes this process 
in a non-perturbative fashion. 



Here, the frequencies uj a p = [E a — Ep + (e a — ep)] /H involve energy shifts e a € R, which are 
defined below. They describe the coherent modification of the system energies due to the 
presence of the environment analogous to the Lamb shift. More importantly, the incoherent 
dynamics effected by the environment is described by the set of complex rate coefficients 
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(a\Tr E (TV 1 / 2 [|a )(/?o|»Pgas] V^V 



(10) 



Our main task is to calculate these quantities. As a fist step we need to specify the rate operator 
r. It is naturally given in terms of the current density operator j = ngasp/n- of the impinging 
gas particles multiplied by the channel-specific total scattering cross sections a (p, a) , 



"-gas ^ n Q <g> —a (p, a) , 



(11) 



with projectors U a = \a)(a\. 

For normalized and separable system-gas states the expectation value of V thus yields the 
infinitesimal total collision probability — provided the motional state of the gas particle is a wave 
packet heading towards the origin. However, since (llip depends only on the modulus of the 
velocity |pj /m it will yield a finite collision probability even if the particle is heading away from 
the origin. Hence, for ([7]) to make sense either the definition of V should contain in addition a 
projection to the subset of incoming states, or the S-matrix should keep such a non-colliding 
state unaffected. The latter is not guaranteed in general, since even a purely outgoing wave 
packet, located far away form the origin, may get transformed by S since the definition of the 
S-matrix involves a backward time-evolution [9]. 



3.1. A phase space-based calculation 

In the following, we will explicitly disregard the outgoing states by using a non-diagonal 
decomposition of pgas- As discussed in [10] the thermal gas state can be written as a phase 
space integration over projectors onto Gaussian states of minimum uncertainty, 

characterized by an inverse temperature (3 greater than the inverse gas temperature (5. Denoting 
the Maxwell-Boltzmann distribution corresponding to the temperature = — j3~ l by 



A(Po) = (27rm//3) 3 / 2 exp(— /3pg/2m) the state flSJ) is represented as 



Pgas = J dp p. (p ) J ^- IVVopqXVVopJ- (13) 
Inserting (fT3|) into (fTOj) yields the expression 

= /dPbA(Pb)/^m5 A (r 0) p d ) (14) 



with 



(ro,Po) := y"dp(a|( P |Tr 1 /2|ao)|^r p D )(A)KV'rop b |r 1/2 Tt|/3)|p). (15) 



Evidently, the phase space function rn^jf (»*o>Po) gives the contribution of different phase space 

regions to the rate coefficient . 

This permits now to restrict the calculation to incoming wave packets. Since the Tn^i are 
averaged over all available positions in (|14p it is natural to confine this spatial average at fixed 
p Q to a cylinder pointing in the direction of p , whose longitudinal support A Pq vanishes at 
outgoing positions. Its transverse base area is given by an average cross section E p . In order 
to evaluate the phase space function (|15|) . one inserts momentum resolutions of unity between 
the T and T operators and notes the representation [9] 

(a/Kp/lTMIft) = Uf Pi > 6 {E Pfaf - E Piai ) (16) 

in terms of the multi-channel scattering amplitude and total energies E pa = p 2 /2m + E a . In the 
calculation one chooses /? large and eventually takes the limit (3 — > oo, — > of very extended 
wave packets so that fi approaches the original Maxwell-Boltzmann distribution [i. In this limit 
the dependence on A Po and E Po drops out provided one identifies S Po with the geometric mean 
of the total cross sections of the involved channels. One obtains [3] 

Kf° = xZ P °^r J d P ^ (Po) Ua (P,P ) fko (P'Po) 6 + E a - E a )j 

(17) 

with the Kronecker-like factor which is equal to one if E a — E aQ = Ep — Ep Q and zero 

otherwise. Moreover, the energy shifts are determined the real parts of the forward scattering 
amplitude, 

e a = -2irh 2 ^ /dpo/i(p )Re[/ aa (p ,p )]. (18) 
m J 

For the special case of factorizing interactions, Hj n ^ = A®B^, and for times large compared to all 
system time scales this result can be obtained rigorously [11] , by means of a "low density limit" 
scaling method [12 [5]- As discussed in [3], limiting forms of ([9]) display the expected dynamics. 
On the diagonal, it reduces to a rate equation, where the total inelastic cross sections determine 
the transition rates, while in the case of elastic scattering the coherences decay exponentially. 



3.2. An equivalent approach 

It is worth noting that the result (|17p for the rate coefficients M®g can as well be obtained in 
a more direct, while less obvious way, if the diagonal momentum representation of Penv is used 
instead of (|13p . In this case, the improper eigenstates \p) have a meaning beyond their use as 
convenient basis states for the expansion of the incoming state. Rather, each momentum state 
in penv must be viewed as representing the limiting form of a normalized wave packet. It is 
therefore necessarily restricted to a discrete set defined by the corresponding volume-normalized 
states \p). The problem is then that the application of S to these spatially extended states 
leads to the unwanted transformation also of its "outgoing components" . As a consequence, the 
resulting expression for M"2^° involving 



is ill-defined since it contains the square of the (^-functions in (|16p and the normalization volume 



This can be healed by the second option mentioned above, the formal modification of the 
operator S, such that it keeps outgoing components invariant. In the continuous basis of 
improper states the unitarity of S is expressed by the optical theorem, which quantifies the 
diffraction limitation of the scattering probability. Since diffraction cannot be accommodated 
in the discrete basis, the requirement of probability current conservation must be incorporated 
additionally in any consistent modification of S, which serves to define a normalized transition 
matrix for the "scattering" into the different "leads" characterized by the discrete \p). The 
point is now that this normalization condition can again be viewed in the continuous basis of 
improper states, and thus provides a simple rule how to form a well-defined expression [10j . Its 
multichannel version is to replace (|19p by 



The appearance of the cross sections in the denominator is here a direct consequence of the 
normalization requirement for the probability current. By using this replacement one obtains 
the result (|17p immediately from (jlOp for any momentum diagonal penv- 

It should be mentioned that the problem of squared delta functions can also be circumvented, 
as discussed in [131 I14j . by appealing to the Fermi golden rule when one writes down the 
dynamic equation. This way one leaves the realm of (time-dependent) scattering theory [9], 
since one needs to speak about the "elapsed time" during the scattering process. Effectively, 
the "interaction Hamiltonian" of the dynamic description is thus identified with the off-shell 
T-matrix, a step which is valid (only) in second-order perturbation theory. It should be noted 
that the form of such a Born approximation to the master equation can differ from the non- 
perturbative result, as demonstrated in the following section. 

4. Collisional decoherence of a Brownian particle 

As a second, nontrivial application of the monitoring approach, let us consider how the motional 
state of a single Brownian particle is affected by a gaseous environment. In analogy to the 
classical case, the corresponding master equation may be called a linear quantum Boltzmann 
equation. Here the term "linear" refers to the fact that the equation for the single Brownian 
particle is linear (the particle is not interacting with itself), and should not be confused with the 
"linearized", i.e., perturbative description of the reduced single particle state a self-interacting 
gas. The resulting equation will describe, on equal footing, both the short-time decoherence 



(2irhf 

n 



(ap 1 1 T 1 a p ) (PoPo I T o I PPi ) 



(19) 




(20) 



behavior, i.e. the rapid "localization" of a spatial superposition state into a mixture, and the 
long-time dissipative behavior, i.e. the gradual thermalization of the Brownian particle. 

Let us denote by M and m the masses of the Brownian and the gas particle, and the reduced 
mass by m*. Since a collision affects only the relative coordinates we denote the single particle 
S-matrix by So = I + iTq. Moreover, it is convenient to denote the relative momentum by 

rel (p,P) := —p-^P- (21) 
m M 

Like above, the first step in obtaining the master equation is to specify the rate operator. 
In classical mechanics, the rate is determined by the modulus of the relative current density 
j re l (p, P) = ng as |rel (p, P)| /m* multiplied by the total scattering cross section a (pj n ) ■ This 
suggests 

T = j rel (p,P)a(rel(p,P)), (22) 

and, indeed, for normalized and separable particle-gas states the expectation value of this 
operator yields the collision rate experienced by the Brownian particle. Like above, the definition 
(f22j) does not yet include the projection to the subspace of incoming particle-gas wave packets. 



4-1- Calculation in the momentum diagonal basis 

Inserting the diagonal momentum representation of the gas state (|8|) into ([7]) we find that the 
equation of motion for the Brownian particle in momentum representation, p (P, P') = (P\p\P'), 
takes the form 

d t p (P, P') = JdQ M in (P, P'; Q) P {P- Q, P' - Q) 

Mo'ut (P) + M Q c i t (P')l p (P, P') . (23) 



1 



In particular, there is no contribution of the second term in (J7]), since the gas density is uniform. 
Equation (|23|) is specified by the complex rate function [2] 

M in (p, p' ; q) = [ d Po p (p ) Vlft + <zk (Pi + q)V\Pi - 9k (Pi - q) 

771* J 

x ^L( pf + g | To | p . + _ g| T t| p/ _ q y (24) 

Here p i := rel ^p , P \ P — Q^j and pj := p^ — Q were introduced as functions of p . Moreover, 
q := rel ^0, P 2 P J . The rate function (f24l) defines also 



Km 



(P) = J dQM in (P + Q,P + Q-Q), (25) 



which will turn out to be the rate in the classical linear Boltzmann equation of a particle with 
momentum P to be scattered by the gas into a different direction or velocity. The remaining 
difficulty is to evaluate the function (f24"|) . Like in the previous section, our use of a diagonal 
representation of the normalized gas state leads to an ill-defined expression involving a square 
of delta-functions and the normalization volume, which can be traced back to the unwanted 
transformation of the outgoing parts of the normalized state, 

^^(Pf + q\T \pi + q){Pt - q\Tl\p f - q)g (q) ■ (26) 



Like above, it can be argued that any consistent modification of the S-matrix, i.e., of To, which 
keeps the outgoing parts invariant must necessarily transform (|26[) into 




f{Pf + Q±,Pi + q±) r{Pf-Qx,Pi-Q±)_ g ^ , (2?) 



y/a (pi + q ± ) \pi + q 1 _\\f<r{pi- q±) \pi - q±\ 



2 



Here, we denote, for given q ^ the parallel contribution of a vector p by p» q = (p ■ q) q/q 
and the perpendicular one by p ±q = p — p» q . These projections manifestly guarantee the 
conservation of energy in the two scattering amplitudes individually. As discussed in [2], this 
treatment leads directly to 



M in 



P,P',Q) = [ dKL (K, P - Q,Q) L* (K, P' - Q, Q) (28) 



with 



xf (rel{K ±Q ,P ±Q ) -Q,rel{K± Q ,P± Q ) + Q\ . (29) 

The integration in (I28p is over the plane = {K € M 3 : K ■ Q = 0} perpendicular to the 
momentum transfer Q. Clearly, the function L contains all the details of the collisional 
interaction with the gas. It involves the elastic scattering amplitude /(Pout'Pin)' ^ ne 
momentum distribution function /x (p) of the gas, and its number density ngas- 

4-2. Operator form of the quantum linear Boltzmann equation 

The particular form (|23p with (|28p and (|29p admits the master equation to be written in a 
representation-independent fashion. One finds the form [2] 



dtp = 7^[H,p] + J dQ J ^dK I^Lq^kP^^ — ^P^q >k ^q,k — ^q K \-q,kP^ 



(30) 



with the Lindblad operators \-q,k defined by 

Lq,k = e iX Q / h L (K, P; Q) . (31) 

Here, the position and the momentum operator of the Brownian particle are denoted by X and 
P, respectively. Reassuringly, this form of the master equation is in agreement with the general 
structure of a translation-invariant and completely positive master equation, as characterized 
by Holevo [T5], see the discussion in [15} H7). Moreover, on the momentum diagonal it reduces 
to the classical linear Boltzmann equation for the momentum distribution function. 

Equation (I30p should be viewed as a full quantum version of the linear Boltzmann equation. 
Although its representation independent operator form (I29p - (|3ip might appear complicated at 
first sight, it has in fact a rather suggestive interpretation when viewed in the momentum 
representation f)23f) . To see this let us define as a possible two-particle momentum trajectory 
(P m ,p m ) — ► (-Pout' Pout) an y combination of initial and final momenta of the Brownian and 
the gas particle, respectively. The allowed two-particle trajectories are those which conserve the 
total momentum end energy. The in-rate (j28[) in the momentum representation f)23[) can now 
be understood as the integral over the scattering amplitudes of all allowed pairs of two-particle 
trajectories, which end at P and P , respectively, and correspond to the common momentum 
exchange Q, cf. Fig[3j The integration measure, i.e., the weight of the pairs of trajectories, is 
simply given by the distribution in the gas, as restricted by the choice of P, P' , and Q, and by 
the requirement of energy and momentum conservation. 
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Figure 3. A pictorial view of the pairs of 
"two-particle trajectories" entering the in-rate f)28|) 
with (|29p . For a given pair of final Brownian 
particle momenta, P and P and for given common 
momentum exchange Q, there are many different 
pairs of (£, £') of initial gas momenta which are 
consistent with the conservation of energy and total 
momentum. The integral in (|28p covers them all and 
weights them according to their distribution in the 
gas. 



4-3. Relation to previous results 

A master equation very similar to the monitoring result (|29p -(j3ip was proposed, on a slightly 
more heuristic basis, in 1995 by Diosi [13J. In the present notation his equation is given by (|3U|) 
and (J2U) provided the functions L (K, P; Q) are replaced by 

Dl0sl V ' V Qml P V v M J 2 M J J \ m 2 ' m 2 J 

(32) 

The main difference with respect to the present result (|3ip is that the arguments of the 
scattering amplitude do not depend on P. This implies that the complex rate function M^ n 
is determined by the differential cross section dcr/dO = |/| 2 , rather the individual complex 
scattering amplitudes of the collision trajectories. This indicates that it does not incorporate the 
non-parallel pairs of trajectories, £ ^ though, a priori, one would expect them to contribute 
to the rate as well. On the diagonal, i.e., for P = P', Diosi's equation coincides with the present 
result (|29p - (|3ip ; however they differ in the "Brownian limit" of small momentum transfers, thus 
predicting different coefficients for the corresponding Fokker-Planck equation. 

Related and more recent works can be found in |X8|, 119] . The most relevant development for 
our case is the theory by Vacchini [20\ [2Tj , which uses the van Hove expression for macroscopic 
scattering to relate the two-point correlation function of the gas with the collision kernel, a 
treatment which is valid in the weak-coupling limit. It turns out that his equation coincides 
with the limiting form of the present result if one replaces the scattering amplitude in (|29p by its 
Born approximation /g. The reason is that /g depends only on the difference of the momenta so 
that the P- and -fC-dependence of the scattering amplitude vanishes in this limit. Incidentally, 
this shows that the non-perturbative form of the master equation cannot be specified by a 
weak coupling calculation alone. In particular, it is not obtained by simply replacing the Born 
approximation of the scattering amplitude with the exact one. 

Another limiting form of the master equation is the case of an infinitely massive Brownian 
particle, m/M — > 0, where it describes no dissipation, but pure spatial decoherence. As one 
expects, the present result (|29l) - (f3TT) reduces in this limit to the proper master equation for 
collisional decoherence [10] . which was recently tested experimentally [22] 123]. 

5. Conclusions 

In summary, I presented a general method of incorporating proper time-dependent scattering 
theory into the dynamic description of open quantum systems. Its derivation is based on the 
theory of generalized and continuous measurements and it yields completely positive master 



equations, which account for the environmental interaction in a non-perturbative fashion. When 
applied to either the case of an immobile, complex object or the case of a Brownian point particle, 
it provides a detailed and realistic account of the dissipation and decoherence effects induced by 
the presence of a gaseous environment. 
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